Fluctuation relations without microreversibility in nonlinear transport 
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In linear transport, the fluctuation-dissipation theorem relates equilibrium current correlations to 
the linear conductance coefficient. For nonlinear transport, there exist fluctuation relations that rely 
on Onsager's principle of microscopic reversibility away from equilibrium. However, both theory and 
experiments have shown deviations from microreversibility in the form of magnetic field asymmetric 
current-voltage relations. We present novel fluctuation relations for nonlinear transport in the 
presence of magnetic fields that relate current correlation functions at any order at equilibrium to 
response coefficients of current cumulants of lower order. We illustrate our results with the example 
of an electrical Mach-Zehnder interferometer. 

PACS numbers: 73.23.-b, 05.40.-a, 72.70.+m 



Introduction - Onsager derived the symmetry of trans- 
port coefficients of irreversible processes using the prin- 
ciple of microscopic reversibility for the fluctuations of 
the equilibrium system 0, 0]. Thus the symmetry of 
transport coefficients in the linear transport regime is di- 
rectly related to the fluctuation-dissipation theorem of 
Einstein, Johnson, Nyquist and Kubo [H, Q, H, [fij . Nat- 
urally, the question arises whether there are fluctuation 
relations beyond the linear response regime. In statistical 
mechanics, fluctuation relations were derived 0, [1| as an 
extension of Onsager's relation to systems far from equi- 
librium. These fluctuation relations make statements on 
the distribution function of observables conjugate to ther- 
modynamic forces for a wide variety of non-equilibrium 
systems 0, Eo| . In electrical transport the variable of 
interest is the transferred charge. The theory is known 

and fluctuation rela- 




as full counting statistics [Hi, [12 
tions for conductors have been discussed in the absence 
IH HI 14, ]j| and the presence of a magnetic field [lj|. 

At equilibrium, in the presence of a magnetic field, 
Onsager reciprocity still holds. However, away from equi- 
librium, the potential landscape inside the conductor is 
neither an odd nor an even function of magnetic field. 
As a consequence, electrical conductors exhibit manifest 
deviations from symmetries based on microreversibility 
and fluctuation relations derived from this principle [161 ] 
are not valid. Surprisingly, and this is a central point 
of our work, we obtain novel fluctuation relations even 
without invoking the principle of microreversibility. Im- 
portantly, the novel fluctuation relations are general and 
independent of a specific model for interactions. 

Full counting statistics and fluctuation theorem - The 
full counting statistics of a conductor with M termi- 
nals is the probability distribution P(Q) that Q = 
(Qi, Q2, ■ ■ ■ ! Qm) charges are transmitted into the reser- 
voirs during the measurement time t. The distribution 
function P(Q) is expressed by the generating function 
F(iA) = ln£ Q P(Q)e iAQ , where A = (A x , A 2 , . . . , A M ) 
are called counting fields. In the long time limit, all irre- 
ducible current cumulants at zero frequency are obtained 



FIG. 1: (color online) Mach-Zehnder interferometer at fill- 
ing factor 2. Only the outer edge state enters the interfer- 
ometer. Coulomb-interactions between neighboring edges -as 
indicated by shading- lead to internal potentials U% and Uo 
in the inner and outer edge. For reversed magnetic field, all 
arrows point in the opposite direction. 



by consecutive derivatives of the generating function, in 
contact a this is <(A/ Q ) fe ) = {-ie) k [d k F/d\ k a ] A=0 /t. 
The magnetic field B perpendicular to the conductor 
and the affinities A = {-j^f, j^f, ■ ■ ■ , f^fO are exter- 
nally controlled. Here eV a is the potential at terminal a 
and T the temperature, assumed to be equal and non- 
zero in all terminals. 



The fluctuation relation for the full counting statistics 
gives a simple relation for the probability that Q or - 
at reversed mag netic field — Q charges are transmitted. 
Derivations [12, EH 14, EH, Ell rely fundamentally on 
microscopic reversibility: a process from terminal a to 
f3 has the same probability as the reversed process, from 
terminal 8 to a at inversed magnetic field. Refs. 
13, 14, [la, [ldl assume that this is valid also far from 



2 



equilibrium and find 



the generating function in the scattering approach is 11 1 



P(Q,B) 
F±(iA,A) 



e AQ P(-Q,-B) 
±F±(—iA - A, A) 



(1) 
(2) 



F(iA) = - I dEtx 



Eq. @ is the Fourier transform of Eq. |T]) and deter- 
mines the symmetry of the generating function. For 
convenience, the (anti-)symmetrized generating function 
F±(iA) — F(iA, B) ± F(iA, ~B) is used, and the nota- 
tion F±(iA, A) emphasizes that the generating function 
depends explicitly on the affinities A. The derivation 
[II H, 0, EH of Eq. J2]) relies fundamentally on mi- 
croscopic reversibility: a process from terminal a to 
has the same probability as the reversed process, from 
terminal to a at inversed magnetic field (see appendix 
A). 

The Mach-Zehnder interferometer (MZI) - As an in- 
structive example, we present a MZI (see Fig. 1) and 
show how interaction (screening) effects lead to devia- 
tions from reversibility. It is a four terminal conduc- 
tor with two quantum point contacts (QPC) acting as 
beam splitters as shown in Fig. [1] The two interferome- 
ter arms enclose a magnetic flux $, such that interference 
arises due to the Aharonov-Bohm effect. In experiments 
17, 18, [3, [5(1, the MZI is realized using edge states 



in the quantum Hall regime, and it is often operated at 
filling factor 2. Only carriers in the outer edge enter 
the interferometer and are able to interfere. Here, the 
inner edge state moves in vicinity in both interferome- 
ter arms and carries current from terminal 1 to 3 and 
from 2 to 4 Although a four-terminal conductor, 

the MZI is characterized by only a single transmission 
probability T31 , due to the separation of left- and right 
movers. T31 is the probability for a particle to be trans- 
mitted in the outer edge state from terminal 1 to 3. In 
the linear transport regime reciprocity means that 21 1 
Tsi(+B) = Tis(—B). We next demonstrate that already 
Hartree interactions lead to a violation of Eqs. |T]) and 
©■ 

Breakdown - Interactions can lead to magnetic field 
asymmetry in nonlinear transport, as was shown theo- 
retically [2S 




as well as experimentally |25|, [2 



1271 . I28t |29j|: Every particle is moving in a local poten- 
tial U(f) generated by all the other particles. The in- 
ternal potential has to be determined self-consistently 
and depends on all potentials Vj applied in the exter- 
nal contacts, U(r) — U(r, {V^}). The scattering ma- 
trix depends on the energy of the particle and is a func- 
tional of the internal potential, S = S(E, U(rj). Indeed 
the functional dependence of the scattering matrix is re- 
quired by gauge invariance: The generating function has 
to be invariant under a shift of all potentials by the same 
amount, V 7 — > V 7 + Uq- This condition can be expressed 
as jjj^ — 0. For long times and neglecting interactions, 



m(l - fK) 



(3) 



Here, K = (1 — A^^AS) is composed of the scattering 
matrix S, the unit matrix 1 and the matrix A introduc- 
ing the counting fields, A = diag(e~ lAl , e _iA2 , . . . ,e~ iAikr ). 
The diagonal matrix / contains the Fermi-functions of 
the different terminals with / = diag(/i, /2, . . . , /m)- 
With this we can show that gauge invariance requires 



^ dK dK n 



(4) 



Here we used that the derivative of the Fermi-functions 
with respect to Uq can be expressed as an energy- 
derivative and that dK/dUo — X^-y dK/dVy. Therefore, 
the scattering matrix entering the expression for K de- 
pends not only on the energy E of the carriers but also via 
the internal potential landscape on the external voltages 
S = S(E, {V 7 }). As mentioned above, the local internal 
potential has to be determined self-consistently and it is 
not necessarily an even function of magnetic field [221 ] . 
As a consequence, for nonlinear transport, the scattering 
matrix is not reversible, S a j3{B, {V^,}) ^ Sp a {—B, {V-y}). 
This implies immediately the breakdown of the fluctua- 
tion theorem ([1]) and ([2]) for nonlinear transport, since 
any derivation is based upon reciprocity. 

The lack of reversibility can be shown explicitly for 
the Mach-Zehnder interferometer. Coulomb-interactions 
between the two edge states moving through the inter- 
ferometer lead to internal potentials U , Ui in the outer 
and inner edge. In this respect the inner edge acts as a 
gate on the outer edge. For the interference, this gives 
rise to an additional phase difference (p(B) = eAU T/h 
between the two interferometer arms. Here, r is the time 
an electron needs to traverse the interferometer, and AU 
is the difference of the internal potential in the outer edge 
between upper and lower arm. 

It is easy to see that the internal screening potential 
U Q is not an even function of magnetic field: For positive 
magnetic field as shown in Fig.[U only processes from left 
to right contribute, and the internal potential will depend 
on the reflection Ra = 1 — Ta of the left beam splitter 
and on the voltages V\ and Vi- For inversed magnetic 
field, processes from right to left are important, which 
depend on Rb = 1 — Tb and voltages V3 and V4. To be 
explicit, we determine the local internal potential self- 
consistently within a Hartree approximation 0,111 El. 
The average charges q and qi in the edges of the upper 
interferometer arm are on the one hand expressed as the 
difference between injected and screened charge, and are 
on the other hand determined by Coulomb interaction. 
For positive magnetic field, this determines the internal 
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FIG. 2: (color online) A graphical representation of the rela- 
tions between response coefficients, compare to Eq. (|17[) . The 
number of circles stands for the order of the cumulant, and the 
number of vertical lines mean the order of the derivative with 
respect to voltage as well as the power of the factor (fcsT) to 
which the response coefficient has to be multiplied. The dif- 
ferent heavy lines (different colors) represent derivatives with 
respect to quantities of different terminals. The summations 
go over all possible permutations, where the number of per- 
mutations is given by m with k the number of circles and I 
the number of vertical lines. Higher order correlations follow 
the same rules and can easily be constructed. 



potentials U , Ui in the upper arm through 



e 2 D(V 1 - Ui) = C(Ui - U ) 



(5) 



e 2 D(R A V 1 + T A V 2 - U Q ) = C(U D - Ui). (6) 



Here, G is the geometric capacitance between the two 
edges, and D is the density of states of an edge state. 
Similar equations hold for the lower interferometer arm 
and for reversed magnetic field. To first order in exter- 
nal voltage, the potential difference AU = ^ a u a V a 
is determined by the characteristic potentials u a = 
[dAU /dV a ] eq .. We find u 3 (B) = Ul {-B) = and 
Ul (B) = R A - e 2 DT A /(2C + e 2 D), u 3 (-B) = R B - 
e 2 DT B /{2C + e 2 D), as well as u 2 (B) = -ui{B) and 
Ui(-B) = -u 3 (-B). 

Using the characteristic potentials, the self-consistent 
transmission probability T 31 = T 3 i(+B,Vi — V 2 ) for a 
particle in the interfering edge to transmit from terminal 
1 to 3 for positive magnetic field is T 3 \ = RaRb+TaTb — 
2^R a RbTaTb cos($-</?) with(^(+B) = eui{+B)r{Vx- 
V 2 )/h. For Ti 3 = Xi 3 (— B, V 3 — V4) at negative magnetic 
field, the additional phase is ip(—B) — eu 3 {—B)r(V 3 — 
V±)/h. The lack of reversibility out of equilibrium is evi- 
dent: 



T 31 (+B, Vi - V 2 ) ± T 13 (-B, V 3 - V 4 ). 



(7) 



This means that the the fluctuation relation ([2]) is strictly 
speaking valid only at equilibrium but has corrections 
for finite voltages. In general, taking into account in- 
teractions beyond the Hartree-level will not reestablish 
reversibility. 

Fluctuation relations for correlation functions 
Including interactions, the fluctuation relation ([2|) for the 
counting statistics is not valid anymore, as shown above 
explicitly within a Hartree-model. Nevertheless, we can 



derive fluctuation relations for current correlation func- 
tions. We emphasize that the following section is general, 
no specific model for interactions is needed. It is useful 
to expand the first few cumulants for eV <C kgT, 



la 

S a (3 



Pi 

7 



0{V 3 ) (8) 

(9) 
(10) 



Up to second order in voltage, the mean current I a in 
terminal a is determined by the linear and nonlinear 
conductance coefficients, and G^ 2 ^. The zero- 

frequency current correlations S a p = (AI a AIp) contain 
equilibrium Nyquist noise and the noise susceptibil- 
ity [3l[ which includes the emergent shot noise. Of 
the third cumulant G Q( g 7 = (AI a AIpAIy) , only the equi- 
librium value C^g is used in the following. All response 
coefficients are obtained from the generating function, 
e.g. G^ 2 ^ 7 = — ield 3, F / dXadVpdV^v / 1, where the index 
means setting A and A to zero. 

(Anti-)symmetrizing the above definitions, both the 
fluctuation-dissipation theorem (for +), and the Onsager- 
Casimir relations (for — ), can be formulated concisely as 



o(0) _ h T(r {l) , r {l) \_ ,q(0) 



(11) 



The next order fluctuation relation connects the third 
cumulant at equilibrium which is odd in magnetic field 
with combinations of the noise susceptibility and nonlin- 
ear conductance coefficients, 



(i) 

a/3 

(2) , M2) 



^Li,f3± + Sfa t a± (12) 



-k B T(G { a 'p 1± + Gp t ' a7± + G^ a p ± ) j - TG^ J 7 ± 

These universal fluctuation relations can be extended to 
any order: A current correlation function at equilibrium 
is expressed by combinations of successive response coeffi- 
cients of lower order current cumulants. They are graph- 
ically represented in figure [5J The first two lines of the 
figure correspond to Eqs. (fTTj) and (fT2f . higher order re- 
lations can easily be constructed. 

The derivation of the fluctuation relations is based on 
the following properties of the generating function: 

F±(-A,A) = F±(0,A) = (13) 
F±(iA,0) = ±F±(-iA,0) (14) 

The first equation defines a special symmetry point at 
iA = —A for which the generating function vanishes, just 
as for A = which originates from probability conserva- 
tion. To demonstrate it, for the case of a non- interacting 
system, we ex pan d Eq. ([3]) in terms of multi-particle scat- 
tering events llj and use the detailed balance for Fermi 



,(0) 
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functions / Q (l - fp) = exp(A a - Ap)fp(l - f a ). Nei- 
ther magnetic field symmetry nor microreversibility are 
needed. For a system with electron-electron interactions 
we start from the definition of the generating function 
F(iA) = ln(V lA( ^e lA Q°^ . Here, Q and Q t denote 

the charge operators at time and time t, and the ex- 
pectation value is taken with respect to the initial state, 
described by a grand-canonical density matrix. At time 
the conductor is decoupled from the reservoirs, and the 
initial Hamiltonian Hq commutes with the charge Qo- 
To derive Eq. (fT3|) . we use that the total energy in the 
system " conductor+reservoirs" is conserved at all times. 
For a detailed discussion see appendix B. We emphasize 
that the identity Eq. (|13p is valid without microreversibil- 
ity. In terms of distribution functions, Eq. (|13|) defines a 
global detailed balance relation, 



Q 



Q 



AQ 



(15) 



valid even if Eq. |TJ) is not true. The second equation, 
Eq. (TT4")) represents the fluctuation relation © at A = 
and is a consequence of microscopic reversibility at equi- 
librium. It follows that even (odd) cumulants at equilib- 
rium are even (odd) in magnetic field as expressed by: 



((Ai a ) k )7 = ±(-i) k ((Ai a ) k y_ 



(16) 



Both functions F±(iA,A) and F±(— iA — A, A) can be 
expanded as Taylor series around A = and A = 0. This 
defines general relations for specific Taylor coefficients, 



d"_F ± (-iA-A,A) 



9A* 



1=0 



(?)«' 



d k F±(tA,A) 



which van- 



ish identically due to Eq. (fl3|) . The last term in the sum 
represents the fc'th derivative of the generating function 
with respect to the counting fields, which is the fc'th cu- 
mulant at equilibrium. Solving the above equation for 
this last term leads to 



<(A/ Q ) fc ) 



k-l 

E 

i=i 



(~k B T) 



k-l 



d k - l ((Ai a y) 



k-l 



eq. 



(17) 

This equation relates a correlation function at equilib- 
rium to a linear combination of response coefficients of 
lower order correlations. Together with Eq. (fT6|) -which 
determines the magnetic field symmetry- it defines new 
fluctuation relations for nonlinear transport. In Fig. [2j 
they are schematically represented and extended to the 
general case of a multi-terminal conductor. 

For the MZI, the fluctuation relation Eq. (fl~2|) can 
be explicitly verified within Hartree. We are concerned 
with temperatures and voltages low compared to the first 
plasma mode of an interferometer arm. Due to the sep- 
aration of left- and right-movers, several response coef- 
ficients vanish, in particular the nonlinear conductance 



^3 3i± as weu as the noise susceptibility Sjj^xi- Also 



(i) 



the third cumulant at equilibrium, C'^^ is zero, be- 
cause the scattering matrix is energy-independent for 
equal length of the interferometer arms. But the co- 
efficients G^3 3 and 3 are finite for — B due to the 
internal potential, and vanish for B. Similar argu- 
ments hold for response coefficients with 1^3. Using 
dg/dU = (4e 3 T/h 2 )^R B T B R A T A sin $, Eq. (H2) simpli- 
fies for the MZI to 



2S ( 3 \ } 3± = k B TG^ 3± = ±k B Tu 3 (-B)dg/dU (18) 



-,(2 



25 3ili± = k B TG^[ 1± = k B Tu l {B)dg/dU. (19) 

The fluctuation relation © which does not account for 
magnetic field asymmetry in screening effects, would re- 
quire that the anti-symmetrized part (— ) of the above 
equations is identically zero [l6j]. Measuring a nonlin- 

(2) 

ear conductance coefficient G\ 33 or noise susceptibility 

£3^ 3 which is asymmetric in magnetic field proofs Eq. ([2]) 
wrong. The fluctuation relations Eqs. (|18p and (|19|) are 
linear in temperature, periodic with the magnetic flux <f> 
and depend on the reflection of the beam splitters; they 
can be experimentally verified. 

Conclusion - We have shown that electron-electron in- 
teractions lead to a breakdown of the usual fluctuation 
relations for the full counting statistics in the presence of 
a magnetic field. The reason is, that interactions can in- 
duce effective deviations from reversibility of scattering 
processes out of equilibrium. Instead fluctuation rela- 
tions can be derived which relate correlation functions 
at equilibrium to response coefficients of correlations of 
lower order. These fluctuation relations are valid even in 
presence of magnetic field asymmetry. 
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by MaNEP and the Swiss NSF. 

Appendix A: Derivation of the fluctuation relation - 
The fluctuation relation, Eq. (p}, which holds in the 
presence of microreversibility, was derived in various sys- 
tems [H, [lH 14, [IH, here we demonstrate it for a 



-,(2) 



non-interacting system described by a scattering ma- 
trix. Most directly, the symmetry relation Eq. ([2]) 
in its non-symmetrized form is derived, F(iA, B) = 
F(—iA — A,—B). The generating function Eq. |3]) is 
F(iA) = tj ^H(iA) with the integrand 



H(iA) =mdet 1 + / (VS^XS - 1 



(20) 



The determinant can be expanded in terms of multi- 
particle scattering events [ll|, [HI , 



H(iA,B)=ln £ |det(s { { b Q } } (i?) 



(21) 



cxp (t ^2 Xa - i e Aq ) n fi n ( x a)- 



7 S{a} 7^M 



5 



For conductors with a single transport mode, {a} denotes 
a set of reservoirs from which particles are injected, and 
{b} is the set of reservoirs into which particles are emit- 
ted. The first sum in Eq. (|21[) runs over all possible sets 
{a} and {b}, and represents all possible, distinct ways of 
scattering a number m of particles, with m ranging from 
to M. The probability that m particles are scattered 
from {a} to {b} is given by | det<S|^| 2 , where the matrix 

sf%) is formed by taking the intersecting matrix elements 
of the columns corresponding to the elements in {a} and 
the rows corresponding to the elements in {&} from the 
scattering matrix S. The products over the occupation 
functions of the different terminals in Eq. (|2ip determine 
the probability that exactly m particles from set {a} are 
injected. The exponent contains all counting fields of set 
{a} and {b}. 

Replacing in Eq. (f2"Tj) the magnetic field B by — B and 
all counting fields i\ a by — i\ a — A a leads to additional 
exponential factors. It is useful to define the set {c} as 
the contacts from which particles are injected but not 
into which they are emitted, the set {d} as those into 
which particles are transmitted but not from which they 
are injected, the set {e} as the contacts from which they 
are injected and into which they are transmitted, and 
the set {/} as the set not at all touched. Then, with 
{a} = {c,e}, {b} = {d,e} and £ {a} = {d,f}, $ {&} = 
{c, /}, the summations and products over elements of 
the different sets can be split up, for example ri 7 e{a} = 

n 7 e{c}n 7G { e } and E ae {b} = J2 a e{d} + E Qe{c} - With 
the help of the detailed balance relation for Fermi func- 



tions, / a (l-/a) 



^ A(l — fa), we nn d for products 



concerning the sets {c} and {d} 

cxp ( e A a - j2 A a ) n a n ( j - a) 

a£{d} a£{c} 76(c) ie{d] 

= n (!-a) n a- 

7 g{c} 7 e{d} 



(22) 



The second and crucial point of the derivation is the 
use of reciprocity S a p{B) = Sp a (—B) which implies 

S{b}(~B) = SS(-B). With this, and rccombining the 
sets of Fermi functions one obtains finally 

H(-A-iA,-B) =ln E | det ( S {a}( B )) 2 ( 23 ) 

{a},{b} 

ex p(* e Xa ~ i e a «) n a n ( i_ a)) 

a£{a\ 7G{fc} 7^{fc} 

which is indeed equal to H(iA, B) in Eq. (|2"Tj) . since the 
sum runs over all possible sets {a} and {6}. 

As discussed in the core of the paper, screening ef- 
fects lead to an internal potential which can be asym- 
metric in magnetic field. Then, away from equilibrium, 
the reciprocity relation is not valid, S a p(E,U(B)) 7^ 



S/3 a (E, U(—B)), and with this the symmetry relation (|2|) 
breaks down. 



Appendix B: Symmetry point of the generating function 
- Eq. (|13[) defines a symmetry point for the generating 
function which is valid for any value of the magnetic field, 
F(-A, B) = F(0, B) = 0. It can be derived without the 
use of microreversibility. For a system, described by a 
scattering matrix, we start with Eq. (f2"Tj) at iA = A 



H(-A) =ln |dct(s {a} 



{b} 



(24) 



{a}, {6} 



ex P (- e a « + e a «) n a n o- - a)- 

a£{a} a£{b} ~f£{a} 7^{a} 

Reformulating the sets of contacts {a} and {&} in terms of 
{c} — {/} as introduced above and inserting the property 
(f2"2"j). the exponential factors will be absorbed into the 
Fermi functions and wc find 



ff(-A) = In E | dct ( 5 w)f II A II (!"A) 



{a},{b} 



7 e{b] 7 ^{b} 



= in E n a n a-A)) =0 (25) 

{6} \ie{b} 7 ^{fc} 

In the first line, the sum over {a} can be performed and 
equals one because of probability conservation. In the 
scattering picture it is thus easy to see that the identity 
F±(0) — F±(—A) = is a consequence of both prob- 
ability conservation and the detailed balance for Fermi 
functions. 

For systems with arbitrary interactions, the generating 
function for charge transfer is defined as an expectation 
value with respect to the initial state [33| 



F(iA) = In 



e -iAQt e iAQo 



(26) 



Here, the vector quantities Qo and Qt denote the charge 
operators in the different terminals at time and time 
t. The initial state is described by the grand-canonical 
density matrix po = e~@ 0+ Q° /Zq, with the partition 
sum Z Q = tr[e^ f3 ^ 0+ACl0 ]. The Hamiltonian H is com- 
posed of the Hamiltonians of all reservoirs and of the 
scatterer, which at time are supposed to be decoupled. 
Importantly, the time evolution operator U(t) contains 
in addition the coupling to the reservoir and interaction 
terms and does not commute with Hq. Inserting the ini- 
tial density matrix into the definition above and using 
the cyclic property of the trace as well as the fact that 



6 



Hq and Qo commute, we obtain 



F(-A) = ln(e A ^e- A ^° 



(27) 



In — tr 

Zq 

In -=-tr 

Zq 

In -=-tr 
Zq 



U\t)e A€ *°U{t)e 



AQ 



e- f3flo U\t)e AClo U(t) 



-PH +AQ p /3H 



U{t)e-P Ao U\t) 



The trace is evaluated by tr[. . .} = J2 n ( n \ ■ ■ ■ \ n )i where 
we choose the eigenbasis of the initial Hamiltonian, 
H \n) = e n (0)|n). The state |n) is characterized by a 
configuration of numbers of particles in all reservoirs and 
the scatterer. We introduce the total energy e n (0) of this 
particular configuration. The important point we can 
make is that scattering processes through the conductor 
leave this energy invariant. To proceed, we consider a 
matrix element of the last four operators in the last line 
of Eqs. (|27p. insert 1 = ^ fc \k)(k\ in the middle and find 

(m\e^°tj(t)e-^°U^(t)\n) = (28) 
= J2( m \e' 3Ila \k(t))e-^ \k(t)\m) 

k 

= J2HHt))(Ht)\n) = s nm . 

k 

For the last step we used that Ho\k(t)) = ek(t)\k(t)). 
This is justified, since the time evolution of state \k) sim- 
ply means a change of the numbers of particles in each 
reservoir. Assuming that scattering processes are instan- 
taneous, it will at any time form a basis for Hq. Then, 
using that the total energy is conserved, ek{t) — £fc(0), 
the matrix element ([28]) is diagonal. This result leads to 
F(-A) = lnir[e-^ + A Q°]/Z = as in Eq. (13]). We 
emphasize again that the identity Eq. (|13| is valid for 
arbitrary electron-electron interactions and without the 
use of microreversibility. Special care should be taken of 
the case when (i) the problem is time-dependent, (ii) the 
temperature is not equal in all reservoirs, and (iii) a bath 
allows energy exchange, e.g. via electron-phonon interac- 
tions. Then, we would have to consider energy currents 
as well and introduce additional counting fields that ac- 
count for the transferred energy. In this case, a similar 
relation can be derived. 
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